function map = GenerateColorScaleJunctions()
    %this one had too many colors
%     T = [0, 20, 117; 8,85,255; 0,237,255; 255,224,0; 255,0,0; 131, 3, 255] / 255;
%     x = [0, 0.025, 0.25, 0.8, 1, 1.5];

    %black and white.
%     T = [0,0,0; 0.2,0.2,0.2; 1,1,1];
%     x = [0, 0.025, 2.5];

%     T = [0, 20, 117; 8,85,255; 0,237,255; 255,224,0] / 255;
%     x = [0, 0.2, 0.4, 0.6];


    T = [0, 20, 117; 8,85,255; 0,237,255; 255,224,0; 255, 0, 0] / 255;
    x = [0, 0.09, 0.2, 0.5, 1];
    
    
    map = interp1(x, T, linspace(0, 1, 3000));
end